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Abstract 

We investigate the possibilities offered by the long baseline experiments T2K, Noua, 
LBNE and LAGUNA for the evaluation of the neutrino mass hierarchy and CP violat- 
ing phase 6cp- We consider a neutrino and antineutrino energy in the interval [0.5,12] 
GeV. It is found that the clearest possible distinction between the two hierarchy sig- 
natures is provided by LAGUNA for an (anti)neutrino energy E^p^^ ~ 0.5 GeV in the 
(Ufj^ —7- i>fj_) — )• disappearance channel. For LBNE at Ep ~ 1 GeV the — t- z^^ 
channel may also provide a distinction, although not so clear, and for Noi/a this may 
be even less clear. These results are essentially the same for the ^23 first and second 
octant solutions. As for 5c p determination, LAGUNA also offers the best chances at 
Ey ~ 0.5 GeV in the — )• channel with 023 in either octant. Regarding nonstan- 
dard interactions, the best possibility for their investigation resides in Noz^a whose 
source-far detector distance can become a magic baseline if the hierarchy is normal 
for channels — )• f^, — )• Df., — )• with energy (8.9-9.1) GeV, (11.6-12) GeV, 
(8.8-12) GeV. If the hierarchy is inverse, magic baselines for Noi^a occur for — )• t'e, 
v^, and with energy (11.8-12) GeV, (8.9-12.3) GeV, (8.8-8.9) GeV. For 

LAGUNA a magic baseline appears for the — )• channel only at (4.7-5.0) GeV in 
inverse hierarchy. We have also investigated a possible complementarity between T2K 
and Noi^a. 



*E-niail: crdas@cftp.ist.utl.pt 
^E-mail: pulido@cftp.ist.utl.pt 



1 Introduction 



The solar and atmospheric neutrino anomahes [I] - |8] have since long ago motivated in- 
terest [9] in long baseline (LBL) neutrino experiments ( |T0] - [H]), aiming at an accurate 
evaluation of some of the neutrino parameters 0. Among these, the values of the mass 
squared differences and mixing angles are by now reasonably well established from the data 
of the solar |1] - [7], atmospheric [8] and reactor experiments [16] - [18]. Still, several issues 
on neutrino parameters are still open, namely the absolute mass, the mass hierarchy and 
the magnitude of 5cPi the CP violating phase. LBL neutrino experiments, either from 
accelerators [TT] - [2], [T^], neutrino factories or beta beams [20], [21] will be essential in 
providing the answer to the two most important questions, namely, whether the hierarchy 
is normal (NH, Amg^ > 0) or inverse (IH, Am|^ < 0) and whether or not CP violation 
occurs in the leptonic sector {sin 5c p 7^ 0). 

The recent measurement of the ^13 mixing angle [16], [18] reduces the parameter de- 
generacy [22] - [25] inherent in the LBL three neutrino analysis, so that the above goals 
have become within reach of the forthcoming experiments. Besides the existing T2K [12] 
and Noz/a [13], the latter expected to start taking data soon, two other major projects are 
being considered: LBNE [19] and LAGUNA/LBNO [26]. The first, whose proposal is under 
current evaluation, aims at sending a neutrino beam from Fermilab to Homestake (~ 1290 
km) and the second favours a baseline from CERN to the Pyhasalmi mine in Finland (~ 
2290 km) with a proposal expected to be submitted by 2014 [26] . 

The present generation of LBL neutrino experiments [TT] - [H] relies on high energy 
neutrino beam facilities where charged pions and kaons are allowed to decay in flight. The 
neutrino (antineutrino) beam composition consists mostly of muon neutrinos (antineutrinos) 
with a small mixture of electron neutrinos (antineutrinos) and a small component of tau 
neutrinos if the beam energy is high enough. An alternative to obtain neutrino beams for 
LBL experiments has been proposed some time ago based on muon storage rings [2H], [2S] 
whereby muons produced from meson decays are stored, accelerated and allowed to decay in 
the long straight sections of a storage ring [30]. These sections are aligned in the direction 
of the neutrino detectors which can be located several thousands of kilometers away. Since 
muons and antimuons decay to almost 100% according to 

the precise content of the neutrino and antineutrino beam can be accurately determined. 

There have been many studies of CP symmetries and mass hierarchy via neutrino os- 
cillations for present and future facilities (see e.g. [31] - [3S]) and the prospects for their 
determination have considerably improved after the measurement of ^13 [TB], [TS]. In the 
present work we take advantage of this fact and explore the possibility for extracting Sep 
and the mass hierarchy from the T2K [T2j, Noz/a [13], LBNE [I9] and LAGUNA/LBNO [26] 

^For a recent review on LBL neutrino experiments see |15j . 
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LBL neutrino experiments. As shall be seen, by conveniently choosing a neutrino energy 
range and source-detector distance, the parameter degeneracy involving the mass hierarchy 
and 5c p can be overcome with the data from a single experiment. In fact to this end one de- 
tector located sufficiently far away from the source may be enough, as already noted in [37] . 
Of course the possible redundancy provided by one or two extra detectors will always be 
welcome for the verification and accuracy of the results [§. 

This paper is organized as follows: we start in section 2 with the derivation of the general 
expression for the oscillation probability in matter with constant density without resorting 
to any approximations. In section 3 we start by comparing the results of our numerical 
calculation with the leading term approximations used in the literature [39], [IQ], [H]. From 
the oscillation probability expression we obtain the biprobability plots [39] for a few relevant 
pairs of oscillation channels in which most of our analysis will be based. Thus in subsection 
3.1 we get such plots for fixed distance, varying energy and several values of 5cp, in the 
cases of normal and inverse hierarchies. It will be seen that within preferred energy ranges 
the relative ordering of the iso-^cp curves remains invariant as the energy changes, while 
for other energy regions at the same distance they successively intersect each other. In 
these cases, for a small change in energy any two neighbouring curves may swop their 
relative position. This entails a less accurate, if at all possible, determination of 5cp in the 
corresponding energy range at that distance. In fact the biprobability curves for constant 
energy as a function of 5c p are closed contours which intersect the iso-5cp curves. Since 
the neutrino energy cannot be known with arbitrarily good accuracy, the two parameters 
(beam energy and distance from source) must be conveniently chosen to ensure that the iso- 
5cp curves do not intersect each other within the energy uncertainty range. Moreover the 
regions of sudden variations of the CP phase along the closed contours should be avoided 
so that 5cp could be determined with a reasonable accuracy. In subsection 3.2 we obtain 
the biprobabilities for fixed energy, varying distance and several values of 5c p also for both 
hierarchies. Closed contours are now obtained for constant distances as a function of the CP 
phase and the same rule as before applies for the requirement of its accurate determination. 
Thus far our analysis will only use the first octant solution for the mixing angle and is 
mainly focused on the Noi^a, LBNE and LAGUNA experiments. It turns out that the latter 
experiment appears to be particularly promising as regards its contribution to providing 
definite conclusions for the mass hierarchy and the 5cp range. At the end of section 3.2 
we provide a comparison between the ^23 first and second octant solutions in the three 
experiments. Section 4 is dedicated to a comparison of T2K and Noz/a emphasizing the 
prospective distinction in these experiments between the signatures of mass hierarchies and 
between the signatures of the first and second octant solutions for ^23- Finally in section 5 
we summarize the present work and draw our conclusions. 

§The idea of two detectors placed at different baselines was first proposed in [38]. However, owing to the 
limited knowledge on the mass differences and mixing angles at the time, a solution to the degeneracies was 
far from reachable. 
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2 The oscillation probability 



The starting point for our neutrino oscillation analysis is the standard matter Hamiltonian 
which in the flavour basis reads 



/ 

ff(^) = U 
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where Am^j = mj — mf, 14 = V^GpNe^., is the Fermi constant and N^k is the electron 
density in the k layer. The latter can be expressed in terms of the mass density pk, the 
electron number density and the nucleon mass mjy as 

Nek = ^. (2) 

niN 

For the cases we study we need only to consider propagation within the Earth's crust, for 
which Pk = 3 g cm'^ and Yfc = 1/2. In eq.([l]) U is the usual leptonic mixing matrix PMNS 
relating the neutrino mass basis to the flavour basi^ 

\u^>=U\ui> (3) 

and in a like manner for future purpose one can deflne a unitary transformation relating the 
mass matter eigenstates to the flavour ones 

\ua>=U'\iy-> . (4) 

Given the deflnition of U (eq.(3)), denoting by the mass basis Hamiltonian operating 
between mass eigenstates |z/j> and < z/j| and H^^^ the flavour basis one operating between 
flavour eigenstates \i'a> and < z/^|, we have 

< z/j|if("Vi >=< Jy^\UH'-"'^U^\iy^ >=< vp\H'^f^\v^ > (5) 

hence 

Hif) = UH^"'^UK (6) 

Since H^'^^ and are related by a unitary transformation their eigenvalues are the same, 
therefore denoting by Vf and V respectively the matrices that diagonalize H^-l'^ and H^"^^ 

VjH^f^Vf = V^H^'^W = Hd (7) 

where Hj:, is the diagonalized Hamiltonian. Using (6) one obtains the relation between 
matrices Vf and V, 

Vf = UV (8) 



^We use the Particle Data Group notation 42 for the U matrix. 
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so that UV also diagonalizes the flavour basis Hamiltonian as can be easily checked from 
eqs.(l) and (6). Moreover applying the definition of U (eq.(4)) 

< u^lH^f^W^ >=< u'0)^H^f^U'\u[ > (9) 

which shows that U' also diagonalizes if-^, thus U' = Vf = UV. 

Since the neutrino Hamiltonian (1) satisfies a Schrodinger like equation 

= H^^fl^f^ /3 = e, /i, r) , (10) 

a state produced as flavour a at the origin becomes, after traveling a distance L 

i^^iL) = S^^mO)i^^{0) (11) 
where Si3a{L,0) satisfies the same Schrodinger equation. So in index notation 



{U%i = {U)^,e-'^^''{U*U (12) 



-z / {U'*y,,iH^f\siU')sjdx 
Jo 

where Aj's are the mass matter eigenvalues. The oscillation probability 

P(z/„^z/^,L) = |^^„|2 (13) 

is then evaluated through 

P(z/, Up, L) = (t/'*);3.e*'^^(f/')«(f/')/3je-'^^^(t/'*),, (no sum in a, /3) (14) 

with the following neutrino parameter values |43] . 

sin^di2 = 0.31, sin'^Ois = 2.4 x 10"^ sin^^as = 0.39, 

Amj^ = 7.6 X 10"^eV"^ Am^^ = 2.4 x 10~^eV^. (15) 
For the sake of comparison we will present one example with the second octant solution for 

sin'^923 = 0.61. (16) 
We will use equation (fT4l) throughout the paper as the basis for our calculations. 

3 Noz/a, LBNE, LAGUNA 

The channels to analyse in this section are the ones for which the three experiments (Noi^a, 
LBNE and LAGUNA) are dedicated, namely the muon neutrino disappearance channel 
{u^ — > u^), the inverse golden channel or electron neutrino appearance (z/^ — )■ Ug) and their 
antineutrino counterparts for normal and inverse hierarchy. We denote these channels as /i/i. 
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yue, their counterparts as fifi, fie respectively and accordingly their oscillation probabilities 
as Pfj,fj,, Pfj,e, Pp,p., Pp.e- As mentioned in the introduction our main results are presented in 
terms of biprobability graphs. 

In the current literature the oscillation probabilities for long baselines are usually calcu- 
lated by resorting to approximations up to first or second order in 'small' parameters such 

as m, m, m, m, m 

Ami, ' A^' ~ ' " 

In particular, as it is known nowadays, the approximation e ^ 613 is rather innacurate. 
Here we will use instead the exact numerical expressions from (14) and show in fig.l the 
comparison between our results and those from the approximations used in the literature. 
In the four top panels of fig.l we plot the oscillation probability P^e as a function of distance 
for eight values of 6c p equally spaced from to 360 degrees and in the two bottom panels 
the biprobability plots for the channel pair fie, fie. Here we consider normal hierarchy only 
and a neutrino energy E=2.3 GeV. Notice the large discrepancies in the probabilities and 
the locations of the intersections from panel to panel. These intersections correspond to the 
magic baselines to which we will return in section 3.2. 



3.1 Fixed distance biprobability plots 

In this subsection we just consider the Noi/a [13] and LBNE [12] experiments and study the 
constant 6cp curves as a function of energy for fixed distance. 

Noua is planned to start taking data in 2013 and consists of a 330 ton near detector at 
the Fermilab site and a 14 kiloton far detector. The latter is a liquid scintillator situated 12 
km off-axis at 810 km in the Neutrino Main Injector (NuMI) beam produced at Fermilab. 
It is dedicated mainly to observe z/^ — )■ oscillation along with its antiparticle counterpart. 
The advantage of an off-axis location is that the neutrino energy is nearly independent from 
the parent meson energy, therefore the beam energy spread is much smaller than for an 
on- axis one Also using the NuMI beam, but not yet finally approved, is LBNE with a 
34 kton liquid Argon far detector at 1290 km distance and a neutrino energy in the interval 
[0.5,12] GeV. 

In each of the four panels of fig. 2 we show the biprobability curves for the pair of channels 
fie, fifi for eight values of Sep equally spaced from to 360 degrees with energy running from 
0.5 GeV to 100 GeV. The top two panels refer to a distance of 810 km and the bottom two 
are for 1290 km which are the source-far detector distances for Noz/a and LBNE respectively. 
The left panels are for normal and the right ones are for inverse hierarchy. In each panel we 
plot the biprobability curves for three energies which are ellipses of large eccentricity: for 
810 km we take E=2.3 GeV which is the preferred Noi/a energy (middle ellipses in the top 
panels), E=l GeV and 0.5 GeV (top and bottom ellipses respectively). In order to keep the 
figures as clear as possible, we choose not to continue the iso-^c-p curves for E<1 GeV. 
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The most favourable energy range for the evaluation of 5cp is the one where the iso-5cp 
curves lie the furthest apart since the phase variation is then the smoothest possible along 
each ellipse. Consequently, as it is seen from all four panels of fig.2, a value near E=0.5 GeV 
is the best choice for 5cp evaluation. For 1290 km (LBNE source-detector distance) we take 
E=7 GeV, 1 GeV and 0.5 GeV (top, middle and bottom ellipses in the bottom panels of 
fig.2). Again we omit the iso-^^p curves below 1 GeV and find that the best sensitivity to 
6cp is for 0.5 GeV, although a double degeneracy remains in its determination due to the 
very large eccentricity of the ellipses. In the limit of increasing energies the constant 5c p 
curves merge and eventually coincide at the point (P^e = 0, Pjj^jj^ = 1) whereas they diverge 
from each other for decreasing energy. This is consistent with the fact that the oscillation 
length increases for increasing energy, so for fixed distance the neutrinos oscillate less as 
their energy increases. The same as in fig.2 is done in figs. 3 and 4 for the pairs of channels 
p,e, Jiji and /xe, Jie respectively. In fig. 3 (for /ie, Jijl) and in the top panels, the middle ellipses 
are for E=2.3 GeV, the top and bottom ones for E=l GeV and 0.5 GeV respectively, all for 
810 km distance as in fig.2. In the bottom panels (1290 km) the top, middle and bottom 
ellipses are for E=7 GeV, 1 GeV and 0.5 GeV. Once again the \so-5cp curves are ommited 
below 1 GeV while 0.5 GeV is the most convenient choice as regards 5cp evaluation, if not 
for a double degeneracy. As for the channel pair /^e, Jie (fig.4) the 0.5 GeV contours are 
again the ellipses with the longest major axes. So this energy is also the most convenient 
choice for 5cp evaluation. It will be seen in the next section that an improved situation can 
be found with data from LAGUNA. 

The possibility of distinction between hierarchies appears realistic so long as the neutrino 
energy can be appropriately tuned, as can be inferred from the comparison between the two 
bottom panels of fig. 3 (/le, Jxji channels at the LBNE source-detector distance). In fact for 
E=l GeV, which corresponds to the middle ellipse, it is seen that Pp,p, ~ 0.48 for inverse 
hierarchy whereas for normal hierarchy 0.62 < P^^ < 0.66. Not so clear a distinction can be 
provided at the Noz/a distance for the same energy, as one can conclude from the comparison 
between the two top panels in figs. 2 and 3. 

In a like manner as figs. 2, 3, 4, fig.5 shows the biprobability curves for the /i/i, JiJi channel 
pair in normal hierarchy for 0.5 GeV, 1 GeV and 2.3 GeV at 810 km. In this case the iso- 
5cp curves are nearly superimposed, so that an unrealistically large experimental accuracy 
would be needed in order to obtain information on 5c p- The same situation occurs for 
inverse hierarchy for the same energies and distance and for normal and inverse hierarchies 
at 1290 km with 0.5 GeV, 1 GeV and 7 GeV. Hence the JiJi channel combination is of 
little use, if any. 

3.2 Fixed energy biprobability plots 

In this subsection we extend our analysis to include the LAGUNA experiment and study the 
constant 5c p curves as a function of distance for fixed neutrino energy from the production 
point up to 2290 km. This is the CERN-Pyhasalmi mine distance where the LAGUNA 
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detector is planned to be installed observing neutrinos produced from CERN. We consider 
two neutrino energies: E=2.3 GeV and 0.5 GeV. Our results are plotted in figs. 6, 7 (for 
2.3 GeV) and 8, 9 (for 0.5 GeV). The top panels in fig. 6 contain the biprobability curves 
for the fie, fifi channel pair, the middle ones for fie, /2/2 and the bottom ones for /le, /ie. 
Left hand panels are for normal and right hand ones for inverse hierarchy respectively. The 
iso-6cp curves are shown for eight values of 6cp equally spaced from 0° to 360° as a function 
of distance, all diverging from a common point at zero distance. The constant distance 
contours are for 810 km (Noz/a, dotted lines), 1290 km (LBNE, dashed lines) and 2290 km 
(LAGUNA, full lines), d. Note that the Noi^a and LBNE contours correspond to the ones 
already shown in fig. 2 for /le, /i/i now superimposed on the iso-6cp curves for varying energy 
instead of varying distance as before. In the top right and middle left panels of fig.6 (/ze, /x/x 
channel pair, inverse hierarchy and fie, Jiji channel pair, normal hierachy) all curves intersect 
at one point. This corresponds to one of the channel probabilities being independent of 5c p 
at one particular distance, the so called magic baseline. Such distance can be evaluated by 
expanding the appearance probability 

P,e = \{U')e^e-'^^\U'*),,\' (17) 

and separating the 5c p dependent and independent terms. We obtain 

P^e = 2si2Ci2S23C23Si3Ci3(As5 + Bcs) + 5cp independent terms (18) 

with 

A = sin{\i — \2)L — sin{\i — Aa)^ + sm(A2 — As)/. (19) 
B = C26ii2[l - cos{Xi - X2)L] - cos{Xi - Xs)L + cos(A2 - Xs)L 

where Ci2 = cose^j,... and L is the distance from the source to the detector. In order to 
ensure the independence of P^e on 5c p one must require the vanishing of the quantity in 
brackets in (18). To this end we note that the three arguments of the sines and cosines in 
the expressions for A and B cannot simultaneously vanish for any possible baseline distance 
and moreover such condition would only be a sufficient one for P^g to become independent 
of 5c p- So one must impose 

tan5cp = - -r (20) 
/\ 

for any 5cp, which requires A = B = 0. In this way for E=2.3 GeV one gets for the first 
three magic baselines for the fie channel in inverse hierarchy Lmagic — 1980 km, ~ 3960 km 
and ~ 5940 km (see fig. 7). Hence the common point of all curves in the top right panel 
of fig.6 corresponds to the shortest one at 1980 km. Analogously, for the Jie, fifl channel 
pair in normal hierarchy, the magic baselines occur for fie, the shortest one being located at 

"The preferred neutrino energy for LAGUNA has not been decided yet. Here we consider the same 
energies for aU three experiments so that the corresponding closed contours appear in the same panel. 
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Lmagic — 2008 fcm, and the following ones at approximately 4020 km and 6030 km. They 
can be easily determined in the same way as in the previous case with the replacements 

U' ^ U'* V ^ -V. (21) 

Magic baselines are particularly useful for the investigation of nonstandard interactions [16] , 

HZ!, m- 

Fig.8 shows the biprobability contours for a neutrino energy E=0.5 GeV and the same 
channels as fig. 6 at the Noz/a (dotted), LBNE (dashed) and LAGUNA (full) distances. For 
simplicity the iso-Scp curves are ommited from fig.8. 

From the inspection of figs. 6 and 8 (neutrino energies 2.3 GeV and 0.5 GeV respectively) 
one can infer the prospects for distinguishing between normal and inverse hierarchies. To 
this end one must compare the left and right panels of these figures. Starting with fig. 6 
(middle panels) it is seen that for LBNE alone (bottom ellipse appearing as a line segment) 
normal hierarchy gives 0.012 < P^g < 0.034, P^^ ~ 0.06 while inverse hierarchy gives 
0.042 < P^g ^ 0.078, Pfifi ~ 0.09. Also for LBNE and appearance probabilities only, the 
comparison of the bottom panels if fig.6 shows for normal hierarchy 0.042 < P^g ^ 0.078, 
0.012 < P^g < 0.034 while for inverse hierarchy such probabilities are interchanged. Another 
typical case can be observed from the top panels of fig.6 with the uppermost ellipses which 
correspond to the LAGUNA experiment. Here it is seen that 0.012 < P^g ^ 0.060 and 
P^e ^ 0.020 for normal and inverse hierarchy respectively, hence some overlap exists in the 
two probabilities. As for the channel one has 0.91 < P^^ ^ 0.95 and 0.97 < P^^ < 1 for 
normal and inverse hierarchies, a difference which will be difficult to detect experimentally. 
A similar conclusion may be drawn from the inspection of the remainder of fig.6. Therefore 
the distinction between the two hierarchies appears to be too difficult to trace with a neutrino 
energy 2.3 GeV. 

In contrast, for E=0.5 GeV the distinction between signatures from each hierarchy 
becomes much more plausible especially with the LAGUNA experiment. In fact from fig.8 
and observing the LAGUNA contours, the two top panels tell us that for normal hierarchy 
0-05 < P^e < 0.14, 0.37 < P^^ < 0.40, while for inverse hierarchy 0.02 < P^e < 0.22, 
0.032 < P^^ < 0.073. Similarly from the middle panels in normal hierarchy 0.01 < P^g < 
0.21, 0.36 < Ppfi < 0.50 while in inverse hierarchy 0.02 < P^g < 0.16, 0.03 < P^^ < 0.06. 
Thus it follows that the best possibilities for hierarchy distinction lie on the investigation 
of channels fifi and flfi at an energy near 0.5 GeV in the LAGUNA experiment. In fact in 
all cases a dedicated observation of fig.8 shows that the difference between the probabilities 
for normal and inverse hierarchies is substantial 

P^^iNH) - P^^(/i/) > 0.30 PfifiiNH) - Pf,f,{IH) > 0.32 (22) 

and moreover it is also seen that the other two experiments ( Noi^a and LBNE) cannot 
offer such a possibility. Since this is the most favourable case for hierarchy determination, 
we have also evaluated this probability difference using the second octant solution for 623 
(eq.(16)). The result is 
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P^^iNH) - P^^ilH) > 0.30 PfifiiNH) - Ppp{IH) > 0.29 (23) 

Furthermore if the hierarchy is normal, inspection of the bottom left panel of fig. 8 shows 
that there are realistic prospects for 6c p evaluation. In fact for LAGUNA, whose contour 
appears in this scale similar to a circle, and for the fj,e channel, we have 0.05 < P^e ^ 0.14. 
In most of this interval, namely 0.06 < P^e < 0.13, the two possible values of Pp,e are 
reasonably apart from each other: either Ppg < 0.04 or 0.17 < Ppg 0.21. In this way the 
ambiguity in Sep can possibly be lifted. If otherwise P^e is found to lie close to cither end 
of the interval, namely 0.05 or 0.14, the value of Sep is also unambigously determined, as 
there is only one possibility in each case: Sep ~ 135" or Sqp ~ 315" respectively. So far 
all results obtained are for the 623 first octant solution (eq.(15)). Still for the energy value 
being considered now (E=0.5 GeV) we show in fig. 9 the same as fig.8 for the second octant 
solution (^23 >45°, eq.(16)) whose results are similar. For larger energies the distinction 
between the first and second octant solutions is even less clear as all contours pertaining to 
the three experiments become closer. We will return to this point in the next section. 



810 km 


Channel 


E (GeV) 


P 


/i/i 


8.9 - 9.1 


0.93 


- 0.94 


fte 


11.6 - 12.4 


(1.75 - 


2.0) X 10~3 


/2j2 


8.8 - 12.4 


0.93 


- 0.96 



Table 1: The energy and oscillation probability for a magic baseline in normal hierarchy at 
Nova (810 km): for each oscillation channel and the neutrino energy range the probability 
shown is nearly independent of 5cp- No magic baseline is found to exist for LBNE and 
LAGUNA in the energy range [0.5,12] GeV in normal hierarchy. 



810 km 


2290 km 


Channel 


E (GeV) 


P 


E (GeV) 


P 


lie 


11.8 - 12.4 


(1.1 - 1.2) X 10^3 








8.9 - 12.3 


0.93 - 0.96 








8.8 - 8.9 


0.92 - 0.93 


4.7 - 5.0 


(7.0 - 9.0) X 10-2 



Table 2: The same as table 1 in inverse hierarchy for Nova (810 km) and LAGUNA (2290 
km). No magic baseline is found to exist for LBNE in the energy range [0.5,12] GeV in 
inverse hierarchy. 

Finally given the source-detector distances for Noz/a, LBNE and LAGUNA, we have 
estimated the necessary neutrino energy for each distance to become a magic baseline and 
the corresponding oscillation probabihty. We analysed the energy interval [0.5,12] GeV for 
which all three detectors are designed. Our results are shown in tables 1 and 2. For each 
oscillation channel we show the energy range where the probability is nearly constant with 
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detectors located at 810 km from the neutrino source (Noi/a) and 2290 km (LAGUNA). 
For LBNE no magic baseline is found to exist in this energy range. We note that for the 
jlfl channel both distances can become magic baselines, whereas for /xe with normal and 
p,e with inverse hierarchy none of the distances is suitable. Since at the magic baseline 
the oscillation probability becomes 6c p independent, any significant deviation from the 
oscillation prediction will be the signature of nonstandard interactions [16], |17], |18]. 



The main objective of the T2K experiment [12] is to discover z/g appearance from z/^. The 
collaboration reported their first results in June 2011 j49j which were later improved ^0] 
and also reported evidence for disappearance [51]. The experiment is a long baseline 
off-axis one with a 295 km source-detector distance from Tokai (J-PARC) to Kamioka (Su- 
perKamiokande) and its peak neutrino energy is around 600 MeV. 

Our aim in this section is twofold: to obtain the predictions for the oscillation proba- 
bilities P^e and P^^ at T2K displaying them in terms of biprobability plots and to explore 
the possible complementarity between T2K and Noz/a. Our results are shown in the four 
panels of fig. 10. We consider normal and inverse hierarchies (left and right panels) and first 
and second octant solutions for 6*23 (top and bottom panels). 

We recall that in the previous section the prospective data from one experiment were 
compared at the same energy for different hierarchies and/or different octants. In no case the 
contours from two different experiments were ever directly compared. Thus it was ensured 
that the comparisons were made at the same point in the oscillation phase. Exploring the 
complementarity between T2K and Noz/a, since the two experiments operate at different 
phases, requires a reduction of the probability to the same phase. Given the fact that T2K 
runs with a peak neutrino energy E=0.6 GeV at 295 km distance one must operate Noz/a, 
whose source-detector distance is 810 km, at a neutrino energy value satisfying 



which gives Eatoi/o — 1-65 GeV. In fig.lO the dotted contours refer to T2K (E=0.6 GeV) and 
the full contours to Noz/a (E=1.65 GeV). Hence observing the left and right panels of fig. 10 
on each row, it is seen that 

{Pi,e)Noua > iPt^e)T2K (uormal hierarchy) (P;,e)7Voz.a < iPf,e)T2K (inverse hierarchy) . 
On the other hand observing the top and bottom panels on each column. 



Although these inequalities apply for any value of Sep, the differences in the probabihties 
are so small that the possibility of ever detecting them in this way is slim. 



4 T2K and Noz/a 




(24) 



(PMM)^o.a > {Pf.f.)T2K (^23 < 45°) (P^^)iVo.a < {P^.^.)t2K (^23 > 45°) . 
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5 Summary and conclusions 



We have investigated the prospects for distinguishing normal from inverse neutrino mass 
hierarchies, 8cp determination and first vs. second octant ^23 solutions with T2K, Noi/a, 
LBNE and LAGUNA long baseline experiments. We examined the oscillation channels 
which are possibly relevant for these, namely /le, /i/i and their antineutrino counterparts. 
Owing to the baseline distances involved, the neutrinos are assumed to pass through the 
Earth's mantle only. The starting point for our numerical analysis is the general formula for 
the matter oscillation probabihty derived in section 2. The discrepancies between the results 
from the leading term approximations existing in the literature based on the smallness of 
the mass square differences ratio or the ^13 mixing angle were also examined. 

The main results of our paper are described in section 3.2 and displayed in fig.8. They are 
suggestive of the importance of LAGUNA and its operation at a neutrino energy around 0.5 
GeV. Otherwise we have found that there are also good possibilities to distinguish between 
hierarchies with other experiments provided the neutrino energy can be tuned to reasonable 
accuracy. To this end the oscillation channels which can provide the best information are 
the muon and antimuon disappearance ones, namely z/^ — > z/^ and — )■ z/^. For a neutrino 
energy E ^ 1 GeV the LBNE far detector will be confronted with a probability P^^ ~ 0.48 
for inverse or a probability 0.62 < P^^ < 0.66 for normal hierarchy (sec fig.3, bottom panels), 
whereas for Noi/a these are 0.65 < P^^j < 0.68 (inverse) and 0.55 < P^^ < 0.56 (normal) as 
can also be seen from fig.3 (top panels). On the other hand for a neutrino energy E ~ 0.5 
GeV the possibility of hierarchy differentiation looks much brighter. In fact if the hierarchy 
is inverse, for the /i/i channel in the LAGUNA far detector one expects 0.03 < P^^ < 0.07, 
if it is normal then 0.37 < P^^ < 0.40 (fig.8, top panels). Moreover for the Jijl channel 
and for LAGUNA with inverse hierarchy we obtain 0.03 < P^^ < 0.06 and with normal 
hierarchy 0.36 < P^^ < 0.50 (fig.8, middle panels). The difference between the oscillation 
probabilities P^^ for normal and inverse hierarchy is thus 0.32 or larger and for P^^ it is 
0.30 or larger (see eq.(22)). Hence the /i/i and /2/2 channels at LAGUNA with a neutrino 
energy E ~ 0.5 GeV seem to be the most promising ones to explore in order to find out the 
mass hierarchy. On the other hand the complementarity between T2K and Noz/a, presented 
in fig. 10, does not offer a clear perspective for the mass hierarchy determination. 

For the ^23 solution in the second octant, we have also checked the possible distinction 
between hierarchies. The chances look almost identical as for the first octant solution, 
although slightly disfavoured in the /i/2 channel (see eq.(23)). Distinguishing between first 
and second octant solutions on the other hand looks difficult as can be seen from the 
comparison between figs. 8 and 9 and from fig. 10. 

As for 5cp determination the LAGUNA far detector offers good chances for its feasability 
on the basis of observing the //e and Jie channels. Again the neutrino energy must be tuned 
to E ~ 0.5 GeV. Then if one is able to experimentally distinguish between probabilities 
P^g < 0.04 and 0.17 < P^g ^ 0.21, a relatively narrow interval for 5cp can be determined 
(see fig.8, bottom left panel). Again, for 6*23 in the second octant the prospects are practically 
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the same. 

Finally, regarding nonstandard interactions, since the baseline distances are a priori 
fixed, these will become magic for a conveniently chosen value of the neutrino energy. This 
energy was evaluated in section 3 along with the corresponding oscillation probability. Given 
this probability, any deviation from such a value is a signature of nonstandard interactions. 
In the appropriate neutrino energy interval for the three experiments, namely [0.5,12] GeV, 
the range to search for magic baselines is E^, > 8.8 GeV for Noz/a and (4.7 — 5.0) GeV for 
LAGUNA. No magic baseline apparently exists for LBNE in the above interval. Detailed 
results are shown in tables 1 and 2. 

To conclude, accelerator based long baseline neutrino experiments offer good prospects 
to discover the 6cp phase range, the mass hierarchy and may be nonstandard interactions 
if they really exist. 
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Figure 1: Four top panels (clockwise): P^e o,s a function of distance from to 7500 
km evaluated with our numerical approach, and with the approximate approximations given 
in l4^ , l4^ , 13^ for eight values of 6cp equally spaced from (P to 36(P . In our case (top 
left panel) constant 6cp lines for 22^ , 27 (P , 31^ coincide with those for 13^, 9(P , 4^ ■ 
Two bottom panels: biprobability plots P^e,Pp,e 810 km (left) and 2290 km (right) from 
the neutrino source. Full, dotted, dashed and dot-dashed contours are obtained from our 
approach, and from the approximate expressions in refs. fj^ , l4^ , fS^ - All results are for 
$23 in the first octant (eq.(15)), normal hierarchy and neutrino energy 2.3 GeV. 
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Figure 2: Biprobability plots for the (fie, fifi) channel pair with energy in the range [0.5,100] 
GeV and 623 <45° (eq.(15)) showing the constant 6cp curves (which merge at coordinates 
(0,1) for large energy) and the constant energy contours. Top panels: 810 km distance with 
(from top to bottom) 1 GeV, 2.3 GeV, 0.5 GeV contours. Bottom panels: 1290 km distance 
with (from top to bottom) 1 GeV, 1 GeV. 0.5 GeV contours. Left and right panels: normal 
and inverse hierarchy respectively. Points marked as ■, 0, •, A, A, V, T are for 6cp = CP , 
4^, 9CP, 13^, 18(P, 22^, 27(P, 31^. 
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Figure 4: Biprobability plot for the (fie, fte) channel pair with energy in the range [0.5,100] 
GeV, 5cp <45° ■ Top panels: the contours for 810 km with 2.3 GeV (dotted), 1 GeV (thin 
full) 0.5 GeV (thick full). Bottom panels: the contours for 1290 km with 7 GeV (dotted), 
1 GeV (thin full) 0.5 GeV (thick full). Due to the figure scale, the contour for 7 GeV in 
the bottom panel can hardly be seen. The constant 6cp curves merge for large energy at 
coordinates (0,0) and the values of 5cp are marked as in figs. 2 and 3. 
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Figure 5: Biprobability curves for constant 6cp for the (fifi,fifi) channel pair at 810 km 
source/ detector distance in normal hierarchy with 623 <45° . Their intersections with the 
constant energy contours are marked. Their closeness prevents determination of 6cp- The 
neutrino energy interval is [0.5,12] GeV. 
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Figure 6: Biprobability graphs for neutrino energy 2.3 GeV (623 <45° , eq.(15)): the curves 
for constant 6cp for varying distance and the contours for Nova (dotted), LBNE (thin full) 
and LAG UNA (thick full). Left and right panels are for normal and inverse hierarchy. Top, 
middle and bottom panels: (fie,fifi), (fie,p,Ji) and (fie, fie) channel pairs. Values of 6cp on 
each contour are marked as in fig. 2: □, ■, 0, •, A, A, V, ▼ are for 6cp =CP, 4^ , 9(P , 13^ , 
18(P, 22^, 27 (P , 31^. Merging of 6cp curves occurs for km distance. 
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Figure 7: P^e (inverse hierarchy, 623 <45°, eq.(15))) as a function of distance from to 
6000 km for eight values of 6cp equally spaced from (P to 36(P . Full, dotted, dashed, dot- 
dashed and dot-double dashed lines are for 6cp=CP , 4^, 9(P , 13^, 18(P respectively. 6cp = 
22^, 27(P, 31^ lines coincide with 13^, 9(P , 4^. The neutrino energy is 2.3 GeV. The 
analogue of this graph for normal hierarchy is shown in the top left panel of fig. 1. 
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Figure 8: Biprobability graphs for neutrino energy 0.5 GeV (623 <45° , eq.(15)): the con- 
tours for Noua (dotted), LBNE (dashed) and LAGUNA (full). The values of 6cp on each 
contour are marked as in fig. 2: □, ■, 0, •, A, V, T are for 6cp =CP , 4^ , 9CP , 13^ , 18(P , 
22^, 27 (P , 31^ . Left and right panels: normal and inverse hierarchy. Top, middle and 
bottom panels: (fie,fj.^), (fie,p,p,) and (fie,Jie) channel pairs. The determination of the mass 
hierarchy and 5cp range appear more favourable in this case than for 2.3 GeV (figlS^. 
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Figure 9: The same as fig. 8 for the second octant solution (623 > 45°, eq.(16)). From 
comparison with fig. 8 it is seen that hardly any distinction can he made between the two 
solutions. 
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Figure 10: Comparison between T2K and Nova: left and right panels are for normal 
and inverse hierarchy, top and bottom ones for 623 in the first (eq.(15)) and second octant 
(eq.(16)). Dotted contours are for T2K and full contours for Nova. Values of 6cp are 
marked in each contour and follow the conventions of the previous figures. 
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